A non-abelian spin-liquid in a spin-1 quantum magnet 
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We study a time-reversal invariant non-abelian spin-liquid state in an 5*17(2) symmetric spin 
5* = 1 quantum magnet on a triangular lattice. The spin-liquid is obtained by quantum disordering 
a non-coUinear nematic state. We show that such a spin-liquid cannot be obtained by the standard 
projective construction for spin-liquids. We also study phase transition between the spin-liquid 
and the non-collinear nematic state and show that it cannot be described within Landau- Ginzburg- 
Wilson paradigm. 
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The concepts of 'order parameter' and 'symmetry 
breaking' have been extremely successful for classifying 
various phases of matter and their phase transitions at 
finite temperature. Interestingly, last twenty years have 
provided us many examples where these notions break 
down in the context of zero temperature quantum phases 
and phase transitions and where one needs a totally 
different framework. In particular, the concepts such 
as 'fractionalization' and 'topological order' have been 
shown to be very useful to describe rich physics of vari- 
ous quantum hall systems, frustrated magnets and many 
other strongly correlated systems [if. 

Our current description of topological order draws 
heavily on the projective construction using various slave 
boson/fermion techniques In this paper we present 
an example of a non-abelian quantum spin liquid with 
bosonic spinous where the Schwinger boson construction 
fails to capture even the mean field state. The spin-liquid 
is obtained by quantum disordering a non-collinear ne- 
matic state in a spin S = \ system on a triangular lat- 
tice. Interestingly, the low energy excitations in the spin- 
liquid have non-abelian statistics! To our knowledge, this 
is one of the few known examples of a time-reversal in- 
variant phase of matter that supports non-abelian exci- 
tations [1, Further, the phase transition between the 
nematic phase and the spin-liquid is also exotic and has 
a very large value of the critical exponent "q associated 
with correlations of nematic order parameter. 

A general Hamiltonian describing a spin 5" = 1 quan- 
tum magnet on an isotropic triangular lattice takes the 
form 



(1) 



<tj> 



Additionally, H may have other short-ranged interac- 
tions consistent with SU (2) spin symmetry such as multi- 
ple ring-exchange or second nearest neighbor Heisenberg 
exchange. 



Consider the Hamiltonian H for K > J > 0. The 
ground state of this Hamiltonian in this parameter regime 
has a three-sublattice nematic order where the nematic 
directors on the three sublattices A, B, and C of the tri- 
angular lattice are orthogonal to each other (say along 
x,y and z respectively) As argued by Tsunetsugu 
and Arikawa [5|, such a state may explain many of the 
features^ of the triangular lattice magnet NiGa2Si, in 
particular, the lack of any dipole moment {S), low tem- 
perature specific heat C(T) and finite spin sus- 
ceptibility at T = 0. The directors of the nematic cor- 
respond to 'hard-axes' i.e. the spins on the three sub- 
lattices fluctuate in the plane perpendicular to their re- 
spective directors such that average value {S{r)) — 0. 
Such a state breaks spin-rotation symmetry while pre- 
serving the time-reversal invariance. In terms of the 
the spin operators, the nematic order parameter at site 
r could be described as a rank- two tensor Qi_i^{r) = 
\ {Sf,{r)S^{r) + S^{r)S^{r)) - jS^^. The director at a 
site r is along the eigenvector of Q^j, that corresponds to 
the zero eigenvalue. 

The above ground state doesn't preserve any contin- 
uous subgroup of the original 5*0(3) symmetry of the 
Hamiltonian in eqn. [T] Thus the low-energy fluctua- 
tions around the ground state consist of three goldstone 
modes. Interestingly, the ground state is invariant under 
the discrete subgroup D2 = R%,R^,R% of 50(3) where 
R% (a = x,y,z) corresponds to a global tt rotation of 
all spins about the three orthogonal axes x, y, z. There- 
fore the order-parameter manifold = SO{3)/ D2- This 
identiflcation of the order-parameter manifold allows one 
to characterize the non-trivial topological excitations out 
of the ground state. We recall that in two spatial dimen- 
sions the fundamental group tti^M) of the order param- 
eter manifold Ai is directly related to the combination 
law for the physical point defects 0] . More precisely, the 
defects are classified by the conjugacy classes of the fun- 
damental group. To calculate the fundamental group of 
M = SO{3)/D2, one notes that the lift of L»2 in SU{2) is 
the eight element non-abelian Quaternion group Q. Thus 
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FIG. 1: The ordered state of the three sublattice nematic 
state. The directors on the three sublattices are shown as 
three orthogonal differently colored rods (black, red and blue) . 
As discussed in the text, one may also define the order param- 
eter space in terms of spinous z which reside at the vertices of 
a bigger auxiliary triangular lattice. We take the vertices of 
this auxiliary triangular lattice (marked as 'x') to lie at the 
centroid of every third triangular plaquette of the original 
lattice. 



is homeomorphic to SU{2)/Q. Since SU{2) is simply 
connected while Q is discrete, using the fundamental the- 
orem on the fundamental group, t:i{A4) — Q 7]. As we 
will shortly see that the two-dimensional representation 
of Q in terms of Pauli matrices would be most relevant for 
our purposes. In this representation the five conjugacy 
classes of Q are given by: 



Co = {I}, Co = {-I} 



±iax, Cy = ±iay, Cz = ±icrz 



(2) 



The class Co corresponds to the trivial class i.e. no 
defect. The class Co corresponds to a 360° disclina- 
tion in two of the nematic directors while the class Ca 
{a — X, y, z) corresponds to defects where there is a 180° 
disclination in all but a-axis directors. 

We are interested in constructing a T = spin-liquid 
state obtained by quantum disordering the nematic state. 
Since proliferating topological defects in quantum mag- 
nets often leads to breaking of various symmetries in 
the paramagnet state, one of the simplest ways to ob- 
tain a spin-liquid is to destroy the nematic state with- 
out proliferating any topological defects. To implement 
this, we use an effective lattice model formulated as a 
gauge theory analogous to the formulation of classical 
spin-nematics as a lattice field theory of a Z2 gauge field 
coupled to a vector field Q. 

The order parameter space M = SO{3)/D2 is equiv- 
alent to a set of orthogonal axes ni, n2 and na at each 
vertex of a bigger triangular lattice with the identifica- 



tion Ua = Rria where R is an element of the D2 group. 
The vertices of this new triangular lattice could be taken 
as the centroids of the triangular plaquettes of the orig- 
inal triangular lattice (figure 1). Equivalentally, using 
the identification SO{3)/D2 = SU(2)/Q, it could be de- 
scribed as a quaternion gauge-matter theory whose imag- 
inary time action is: 



S 



QijZj 



A 



^Tr (qij Qjk qki qii) 



Tr (q^j qjk qki) 



Sb 



(3) 



Here z — [z^,z^]^ is a two component spinor with 
the constraint z'^z — 1 which is minimally coupled to 
a quaternion gauge field qij. z's and q's live at the ver- 
tices and the links respectively of a stacked triangular 
lattice. The first term in the action S is the kinetic en- 
ergy term for the spinons z while the second term is the 
kinetic energy term for the gauge fields at a spatial (A) 
or space-time (□) plaquette. The term 5*5 is the Berry 
phase term associated with topological defects. Since 
the topological defects are gapped in all the phases con- 
sidered in this paper, Sb can effectively be ignored for 
our purposes. The original 5*0(3) symmetry of the spin- 
Hamiltonian H is realized as a SU{2) symmetry acting 
on the spinor z. This also indicates how to relate the 
spinor z's to the axes ni,n2 and na. Let us define an 
SU{2) matrix U build from z's: 



U 



(4) 



The 3x3 matrix R build from U 



(5) 



has ni,n2 and ns as its first, second and third 
columns respectively. Thus, all three directors n^'s are 
quadratic in z's and hence the nematic order parameter 
is of fourth degree in z's. 

Let us describe the phase diagram corresponding to 
the action 5* in eqn. [31 First consider S at finite temper- 
ature. Since a continuous symmetry cannot be broken in 
two dimensions with short-range interactions, no symme- 
try broken phases can exist at any non-zero temperature. 
Further, since a pure quaternion gauge theory in two spa- 
tial dimensions is confining at any non-zero temperature, 
in terms of the original spin model only a paramagnet 
phase without any topological order can exist at T > 0. 

Next consider the phase diagram at zero temperature. 
Clearly for t ^ k, the spinons would condense yielding 
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a three sublattice nematic state. As discussed above, 
this phase would have three goldstone modes. The non- 
abehan fluxes through plaquettes would exist as excita- 
tions in this phase and correspond to the non-abelian 
disclinations in eqn. [21 These fluxes have logarithmic in- 
teraction with each other which is mediated by the spin- 
wave magnons. Now consider quantum disordering this 
phase without proliferating these defects. The resulting 
phase would be a paramagnet that would correspond to 
the deconfined phase of a pure quaternion gauge theory. 
The excitations in this phase would correspond to gapped 
spinous, quaternion fluxes and their composites ('dyons') 
which would have interesting fusion and braiding proper- 
ties. Additionally, the ground state of this phase would 
have topological degeneracy on a torus. 

The topological properties of this phase are well un- 
derstood and a detailed discussion of the particles and 
their topological interactions was provided by Propitius 
and Bais in reference Q. The distinct magnetic fluxes 
are in one-to-one correspondence with the topological de- 
fects of the nematic phase, hence they are classifled by 
the conjugacy classes of the quaternion group listed in 
eqn. [2l As in any gauge theory, the electric charges are 
labelled by the irreducible representations of the gauge 
group. In particular, the spinous z-f, transform un- 
der the two-dimensional representation of the quaternion 
group while the flux composites of the form \i(Ta, —i(7a > 
with zero net flux also carry charge and transform un- 
der the one-dimensional representation of the quaternion 
group. Finally, one could have dyonic particles which 
correspond to the bound state of a flux and an electric 
charge. The dyons corresponding to a flux $ are classi- 
fled by the irreducible representations of the centralizer 
corresponding to conjugacy class of $ (one may recall 
that the centralizer of a conjugacy class is the set of ele- 
ments which commute with all elements of that conjugacy 
class). Since the centralizer of the Cq conjugacy class is 
the whole group Q, the bound state of a spinon and a Co 
flux is classified by the irreducible representations of Q. 
On the other hand, the centralizer of the Cx, Cy and Cz 
conjugacy classes is the Z4 group, hence the correspond- 
ing dyons fall into four distinct classes corresponding to 
the four onc-dimcnsional representations of the Z4 group. 

Next consider the ground state degeneracy on a torus. 
Since the deconfined phase is gapped, one could calcu- 
late the degeneracy on a torus by restricting the system 
to a single plaquette with periodic boundary conditions. 
Thus, the number of ground states is proportional to the 
number of inequivalent flux configurations that satisfy 
Y[\j cTijCjkO'kia'ii = 1. A simple counting shows that this 
number is 22. As is well known, the ground state de- 
generacy for a topological phase is equal to the number 
of distinct particle excitations of the system. Thus there 
must be 22 particles in a quaternion gauge-matter theory 
as could be easily verified by direct calculation 

Since spinous transform under a two-dimensional rep- 



resentation of the quaternion group, it's also interesting 
to consider the effect of transportation of a spinon around 
a magnetic flux. For example, when a spinon z carrying 
spin- up, z = [1 0]^ , goes around a flux Cx, it is trans- 
formed to z' = ax z = [0 1] i.e. a spinon carrying down- 
spin. This may seem counterintuitive as it seemingly vi- 
olates global spin-conservation. Since the ground state 
is a spin-singlet, one needs to be careful while discussing 
individual spinons. An instructive thought experiment 
which resolves this paradox is the following. Consider 
putting our system on a torus with a ax vortex thread- 
ing through one of the holes of the torus. Next an up 
spinon, down spinon pair is created out of the vacuum 
and the up-spinon is transported around the vortex keep- 
ing the location of the down spinon fixed. Finally, the up 
spinon is brought back to its original location. The in- 
teresting question is whether the spinon pair remain in 
the spin singlet sector (and hence can be annihilated) at 
the end of this process. 

To answer this, we note that the gauge invari- 
ant wave-function for the pair is given hy \ip) = 
zl.^qm(lrir2---Qrjr'zl./^\Q) whcrc r, r' are the locations 
of up and down spinon respectively and gr„r„_,.i 's are the 
values of the gauge fields on the links connecting the two 
spinons. For concreteness, let us assume that a Cx vor- 
tex threads along the x-axis and is parameterized by the 
gauge choice qrr' = if^x for r = {0,r x) ,r' = r + y Vr 
while Qrr' = 1 otherwise. As the up spinon passes 
through the x-axis, the field operator z-f transforms to z^ 
while the string operator S = qrriQrir2---<lrjr' connect- 
ing the two spinons transforms as 5 — >■ axS = ax since 
in the ground state the only contribution to S comes 
from the vortex gauge field. This implies that at the end 
of the adiabatic transformation, the up-spinon becomes 
a down-spinon and vice-versa since the operator iS acts 
on the down-spinon and transforms it to an up-spinon. 
Hence the spinon pair behave like an EPR spin-singlet 
[lo| and the spin is delocalized while they are separated. 

Having described various possible phases of the action 
S, it is worthwhile to compare our approach to standard 
methods for obtaining spin-liquid states. Our current 
understanding of spin-liquids with bosonic spinons draws 
heavily on the slave-boson formulation. In this approach, 
one represents the spin-S' operator in terms of a two com- 
ponent boson b = [61 62]'^ as S = ^b^crb. The Hilbert 
space of the bosons is projected to the physical subspace 
of a spin by the relation b^b = 25. Clearly, all physical 
operators are invariant under the U{1) local transforma- 
tion, 6(r) — e^'^^'^\r). The physical symmetries of the 
underlying spin-Hamiltonian are realized as transforma- 
tion of spinons b under these symmetries combined with 
a local C/(l) transformation. The spin-liquids correspond 
to states that do not break any physical symmetry and 
their low energy description consists of spinons b coupled 
to a gauge field whose gauge group is always a subgroup 
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of [/(I). 

Could the non-abelian quaternion phase described in 
this paper be described by such a projective construc- 
tion? Interestingly, the answer is no! First and foremost, 
as we mentioned above the standard Schwinger boson 
technique can only describe abelian spin-liquids. Fur- 
ther, the spin-operator cannot be written as an operator 
quadratic in spinous z for the following reason. The ne- 
matic order parameter Q^^y is quartic in z's and since 
it transforms as a spin-two operator, one needs to take 
it's tensor product with itself to construct a gauge in- 
variant physical spin operator S. Explicitly, S is given 
by Sa — —i£°'^'^QbdQdc- Therefore S would be of degree 
eight in z's, which is very different than the usual slave 
boson theories. 

The phase transition between the three sublattice ne- 
matic phase and the non-abelian spin-liquid has also 
many interesting features. This phase transition corre- 
sponds to Higgs transition for the spinous z and as we 
argue, is not describable within Landau- Ginzburg- Wilson 
paradigm. Let us start by analyzing the symmetries of 
the critical action. The symmetry under physical spin 
rotation corresponds to the left multiplication of the ma- 
trix U (defined in the eqn. |4|) by an SU{2) matrix. Re- 
markably, the action is invariant even under the right 
multiplication of U by an SU{2) matrix and thus the 
critical theory has in fact 0(4) ^ SU{2) x SU{2) sym- 
metry! This is because the under unit translation Ta on 
the triangular lattice, the transform as: 



To, -- ni n2, n2 ns, "3 



ni 



(6) 



One finds that the operation of right multiplication of 
U by an SU{2) matrix corresponds to the rotation of ria's 
amongst each other. The generators of these rotations 
Ka (a = 1, 2, 3) satisfy 



(7) 



Since the critical action must have both spin-rotation 
and translational invariance, it would be invariant under 
the transformation U VrUVl where Vr,Vl e SU{2). 
Hence the critical theory is 



S ^ -J d'xTr [(a^C/t) (d^U)] (8) 



- J Sx dr \dfj,Za 



with z'l'z = 1. We emphasize that the critical theory 
can be written in terms of spinous z only because the 
topological defects are suppressed which renders them 
single- valued. We note that the critical theory is very 
similar to that for the phase transition between a spiral 



anti-ferromagnet and a Z2 spin-liquid in a spin S ^ 1/2 
triangular lattice magnet 

The fact that that the nematic order-parameter Q^,y is 
of fourth degree in z has dramatic consequences for the 
critical correlations. For example, the critical exponent 
rj defined by 



{Qtj.u{k,uj)Qf_,iy{k,uj)) 



1 



(a;2 - k^y-v/^ 



(9) 



would have a large value which equals 77 = 3 in a large 
N limit if one generalizes the 0(4) model to an 0{N) 
model. For finite N, one would obtain 77 > 3 whose pre- 
cise numerical value we do not calculate here. This is very 
large compared to the anomalous exponents correspond- 
ing to the order parameter in usual Landau-Ginzburg- 
Wilson theories. 

Apart from the order parameter, the six conserved cur- 
rents associated with the 0(4) symmetry would also have 
power-law correlations at the critical point. Three of 
these, the conserved total spin Stot = J2r '^i''') ^^'"^ 
served microscopically while the other three are the KaS 
defined above which are conserved only the low energy 
effective theory. 

The conserved currents acquire no anomalous dimen- 
sions and hence have scaling dimension d = 2. Therefore 
their correlations at the critical point are given by: 



{Ja{r,T)Ja{O,0)) 



1 



(10) 



where J = {S, K}. Comparing eqn. [Slwith eqn. [TUl 
one notices that the conserved currents have a slower 
decay than that for the order parameter which is rather 
unusual. 

One may ask what interactions one might add in the 
Hamiltonian H so as to destroy the nematic state? One 
simple way to obtain a paramagnet state out of the ne- 
matic state is to consider a stacked triangular lattice 
spin S — 1 system and add an antiferromagnet inter- 
layer interaction J±. When J± ^ J^^,K^^, the intra- 
layer couplings, one would obtain a translationally invari- 
ant paramagnet state corresponding to decoupled spin 
S = 1 chains. On the other hand, when J± ^ J|| , i^n , 
one obtains ordered three-sublattice nematic state. Thus 
one can quantum disorder the nematic state by tuning 
J_l/J||, J±/X||. It would be interesting to explore the 
possibility of realizing stacked copies of the non-abelian 
spin-liquid described in this paper in such a setting. 

In summary, we have described a non-abelian spin- 
liquid in a spin S = 1 quantum magnet on triangular 
lattice which can not be accessed within the standard 
slave boson/fermion projective construction. The non- 
abelian phase has interesting topological features cap- 
tured by the ground state degeneracy on a torus and 
braiding and fusion of its excitations. We also described 
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a non-Landau phase transition between this spin-hquid 
and a non-colhnear nematic state. The nematic correla- 
tions near this phase transition are characterized by large 
anomalous dimension. 
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